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Abstract. We study the symmetry of large deviation functions associated with time- 
integrated currents in Markov pure jump processes. One current known to have 
this symmetry is the fluctuating entropy production and this is the content of the 
fluctuation theorem. Here we obtain a necessary condition in order to have a current 
different from entropy with a symmetric large deviation function. This condition is 
related to degeneracies in the set of increments associated with fundamental cycles 
from Schnakenberg network theory. Moreover we consider 4— states systems where we 
explicitly show that non-cntropic time-integrated currents can be symmetric. We also 
show that these new symmetries, as is the case of the fluctuation theorem, are related 
to time-reversal. However, this becomes apparent only when stochastic trajectories are 
appropriately grouped together. 



1. Introduction 

Large deviation theory [IHl] is the branch of mathematics that deals with exponentially 
decaying probabilities. Therefore, it is the appropriate mathematical theory for 
statistical physics. For example, two basic concepts in large deviation theory are that 
of a rate function (or large deviation function), which gives the rate at which a certain 
probability distribution decays exponentially, and the scaled cumulant generating 
function, which is the Legendre-Fenchel transform of the rate function. In equilibrium 
statistical physics, the microcanonical entropy is a rate function and the canonical free 
energy is the corresponding scaled cumulant generating function. Moreover, it can be 
shown that the maximum entropy and minimum free energy principles are consequences 
of the more general contraction principle, which is central in large deviation theory. 
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Furthermore, large deviation theory also plays a prominent role for systems out 
of equilibrium. A series of relations known as fluctuation relations [5H2D] are the most 
general statements known in nonequilibrium statistical physics, for which a general 
theoretical framework still lacks. Among these relations are the Jarzynski relation [5], 
the Crooks relation [13], the Hatano-Sasa relation [H] and the Gallavotti-Cohen-Evans- 
Morriss (GCEM) fluctuation theorem [5H8]. These statements are about very rare 
events, where a certain fluctuating entropy takes a negative value. Therefore, they 
are appropriately described with the use of large deviation theory. In particular, here 
we focus on the GCEM fluctuation theorem for Markov pure jump processes. This 
relation is written as a symmetry in the large deviation function (or the corresponding 
scaled cumulant generating function) associated with the probability distribution of the 
fluctuating entropy and is also known as the GCEM symmetry. 

More broadly, one can consider a class of functionals of a stochastic path known 
as time-integrated currents. Whereas the entropy, which is a specific time-integrated 
current, has a symmetric large deviation, in general other currents do not display this 
symmetry. The GCEM symmetry has its roots on the fact that the entropy is a very 
special functional given by the logarithm of the weight of the path divided by the weight 
of the time-reversed path. Hence, one question that arises is that if it is possible to find 
other currents with a symmetric large deviation function and what would be the physical 
origin of the symmetry. 

This problem has been recently addressed in [22] where it was shown that for a 
restricted class of Markov pure jump processes a current different from entropy presents 
a symmetric large deviation function. As examples, it was shown in [22] that besides the 
entropy, the height in a surface growth model (see also [23]) and the mechanical work in 
a toy model for a molecular motor display such symmetry. More clearly, Markov jump 
processes are used to describe a large amount of physical processes out of equilibrium and 
time-integrated currents are important physical observables. Therefore, a more complete 
theory about symmetries of large deviation functions associated with time-integrated 
currents might be relevant for the theoretical understanding of nonequilibrium statistical 
physics. 

We pursue this direction in the present paper where we obtain a necessary condition 
in order to have the symmetry in a non-entropic current. This condition is phrased in 
terms of Schnakenberg network theory [2T], where the states of the Markov process 
form a network, with the transition rates representing the edges and the states the 
vertices. We show that in order to have a symmetric current different from entropy, 
the current increments related to cycles in this network have to be degenerate. This 
result comes from an analysis of the characteristic polynomial of a modified Markov 
generator [TTlfTTlfTH]. 

In addition, we consider explicitly a 4— states system with three cycles. We show 
that in this case two other symmetries, different from the GCEM symmetry, can be 
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found (see Fig. [5]). We also demonstrate that in a fully connected 4— states system many 
other symmetries arise. As for the physical origin of the symmetries, following previous 
work [22], we show that these symmetries are also related to time- reversal. However, 
they come from the time-reversal of a group of trajectories, where the grouping depends 
on the specific current under consideration and is related to the degeneracies in the 
increments of the cycles. 

We note that links between Schnakenberg network theory and fluctuation relations 
have been addressed previously in the literature [15l[T7l|2ll[25]. The main differences 
between these works and the present paper are: [T7l[2^ study multidimensional joint 
distribution of currents while here we consider the probability distribution of a single 
non-entropic current; we focus our study on elementary currents (see section 3) while the 
work explained in [TJ] concerns fluctuating currents associated with topological cycles 
(Kalpazidou [26] provided a survey of the interconnection between topological cycles and 
edges). Finally, the effect of the coarse graining procedure on the entropy production 
and its relation to Schnakenberg network theory is analysed in [25] . 

The paper is organized in the following way. In the next section we define time- 
integrated currents and show how their scaled cumulant generating function can be 
obtained from a modified generator. In section 3 we introduce Schnakenberg network 
theory, which is used in section 4 where we obtain the necessary condition. Section 5 
contains an analysis of the 4— states and three cycles system. The new symmetry as a 
result of the time- reversal of a group of paths is discussed in section 6. We conclude in 
section 7. Moreover, in the Appendix A an extension of the Kirchhoff 's law in the large 
deviation regime, which is important for Schnakenberg network theory, is presented and 
the 4— states fully connected system is analyzed in Appendix B 



2. Current probability distribution and modified generator 

Pure jump continuous time Markov processes [27] are defined by the transition rates 
between a pair of states. If we consider a pair of states (x, x'), the probability per unit 
time of going from state x to x' (the transition rate) is denoted by w^^^'- In this paper 
we restrict to a countable set of states represented by F. The probability of being in a 
state X G F at time t follows the master equation, which reads 

^P(x, t) = -A(x)P(x, t) + Y^ Pi^'^ t)w^'^x, (1) 

where X{x) = J2x'=/=x "^x-^x' is the escape rate from state x. This equation is also normally 
written in the matrix form = PtC, where C is the Markov generator. It is defined 

as 



if x' 7^ X 
if x' = X 



(2) 
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The object we study here is the so-called time- integrated current (or just current). 
This is a functional of the stochastic trajectory, which is a sequence of jumps 
x{0),x(ti), . . .x(tM-i),x{tM), taking place at random times tj, within a fixed time 
interval [0,T] , where M is the fluctuating total number of jumps in the trajectory. 
More precisely, the stochastic trajectory starts at state x(0) at time t = and then 
jumps iteratively from state to state x{ti) at time t = ti until it reaches a state 

x{tM), where it stays until at least t = T. Representing a stochastic trajectory by ^[o,t]5 
a time-integrated current is a functional written as 

'/T[^[o,r]] = (3) 

i=0 

where f{x, x') is the increment of the current when the trajectory makes a jump x — )■ x'. 
Furthermore, a current is a functional such that this increment has the property of being 
antisymmetric, i.e., f{x,x') = —f{x',x). 

Since the current is a functional of the stochastic trajectory we can consider a 
probability distribution of currents in the following way. Given a time interval T the 
probability that the current ([3]) takes the value aT is written as 

P Ut = a)=Yl P(^[o,T])5( Jt[X[o,t]] - aT), (4) 

where jx = Jt/T is the time-averaged current and F{X[q^t]) is the weight of the path 
X[o,T] • In most of this paper we do not carry the explicit dependence of the current on 
the stochastic trajectory X^q^t], writing only Jt- 

The fluctuation theorem is a symmetry in the probability distribution of the entropy 
current St, which is defined by the increment f(x,x') = In ^ i.e., 

St[X[o,t]\ = 2^ In • (5) 

We point out that we are considering processes such that if Wx-^x' 7^ then Wx'-^x 7^ 0, 
otherwise entropy cannot be defined. This symmetry, known as the GCEM symmetry, 
is valid in the limit of T — )■ oo and is related to events in which the entropy current 
considerably deviates from its average. Therefore, it is conveniently written in terms of 
the large deviation function /^(a), which is defined by 

/^(a) = lim In P(st = a) , (6) 
where st = St/T. Explicitly, the GCEM symmetry reads 

h{a) ~ h{-a) = -a. (7) 
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Note that we are using the subscript s to denote the large deviation function associated 
with entropy. For a general current of the form ([3]), we denote the large deviation 
function by J (a). 

Instead of the probability distribution of a current we can work with the associated 
generating function. Particularly, the scaled-cumulant generating function related to 
is defined by 

i{z) =lim hn{exp{zJT)), (8) 

where the brackets indicate an average over all stochastic paths. It follows from the 
Varadhan theorem [IHl] that I{z) is the Legendre-Fenchel transform of /(a), that is, 

I{z) = sup„gK{^« - ^(a)}- (9) 

Therefore, the GCEM symmetry can also be written as 

h{z) = i,{-l-z). (10) 

More generally, for any current proportional to ^ in the large deviation regime we 
obtain the symmetry I{z) = I{—E — z). Moreover, it can be shown that I{z) is given 
by the maximum eigenvalue of a modified generator associated with the current Jt [TT] . 
This modified generator is defined as 

, , _j Wa:^x> exp{zf{x,x')) H X ^ x' , , 

I -A(x) ifa; = a;' ' ^ ^ 

Note that this is not a stochastic matrix, but it is still a Perron- Frobenius matrix: it 
has a unique real maximum eigenvalue which gives I{z). 

In this paper we are interested in finding currents following the symmetry I{z) = 
I{—E — z) that are different from entropy in the limit of T — oo (i.e. not proportional to 
St)- a sufficient condition for that is a fully symmetric spectrum of eigenvalues [TTtlTS]. 
In this case, the characteristic polynomial associate with C{z), 

P{z,y) = det{C{z)-yId), (12) 

where Id is the identity matrix, follows the symmetry P{z, y) = P{—E—z, y). Therefore, 
our objective is to find currents different from entropy in the large deviation regime such 
that the characteristic polynomial of their associated modified generators is symmetric. 
As a general result in this direction, in section H] we obtain a necessary condition for this 
symmetry to arise in a non-entropic current. Before going into that, in the next section 
we define elementary currents and introduce Schnakenberg network theory, which are 
important in analyzing the determinant (fT2l) . 
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Figure 1. Network of states with 4— states and three cycles. On the right, the spanning 
tree is in blue and the two chords are in red. The fundamental cycle associated with the 
chord (a, b) is ci = (a, b, d, a), while the one related to the chord (a, c) is C2 = (a, c, d, a). 



3. Elementary currents and Network theory 

We now consider the space of states F as a graph where the vertices are the states and 
the edges represent the transition rates. Therefore, if the transition rate between two 
states is zero there is no edge connecting these states. We denote this graph by G'(r). 

Given a pair of states (x, x'), the elementary fluctuating current from x to x' is 
written as Jt{x,x'). Moreover, an elementary current is such that f{x,x') = 1 (which 
implies f{x',x) = —1 ) and the increment is zero for all other pairs of states in G{r). 
Therefore, the general current ([3]) can be written as 



x,x' 



An important restriction on the number of independent elementary currents is the 
finite time Kirchhoff's law 128], which reads 



^ Jt(x,x') = ±1,0. (14) 

This relation comes from the fact that, during a stochastic trajectory, when the system 
reaches the state x, if there is a subsequent jump, the system will leave x. Defining 
jT{x,x') = Jt{x,x')/T, the above relation can be written as 

J2Mx,x') = 0{^) (15) 

This formula is valid for any trajectory (remember that jrix, x') is a functional of the 
trajectory -^[o,t]) and in the typical regime (T — )■ oo), it gives the usual Kirchhoff's law. 
As we explain in Appendix A it is also valid in the (non typical) large deviation regime. 



which is normally not appreciated in the literature (see [22] for a counter example). 

Let us now introduce Schnakenberg network theory [21] (see also [T7t[2Fj). First 
we introduce the concepts of cycle, fundamental cycle, spanning tree, and chord. A 
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cycle in the network G{r) is a closed path (or loop): this is a a sequence of jumps 
c = [xi,X2, ■■■,Xi], which finishes in the same state it started and does not go trough the 
same state more than one time. Note that, except for cyclic reordering, the order of the 
states is relevant. 

Given a current of the form ([3]), the increment related to a cycle is 

n(c) 

K{c) = Y,fi^^,^^+l), (16) 
i=l 

where n(c) is the number of states in the cycle. Furthermore, the product of rates of 
the cycle c, which we refer to as the rate of the cycle, is given by 

n(c) 
i=l 

The set of all cycles in G{r) with at least three jumps is denoted by O = {c/n(c) > 3} 
(note that K{c) = if c is a transposition, a cycle with two states). Given a cycle 
c = [xi,X2, ...x„(c), Xi] , its reverse [xi,Xn{c), ...X2,Xi] is represented by c and, therefore, 
K{c) = -K{c). 

A spanning tree is a undirected connected set of edges that goes through all the 
vertices in G{T) and has no cycles (see Fig. [T]). After choosing a maximal spanning 
tree, all edges that are not part of it are called chords. A chord is represented by / and 
connects the states xi and x'l (from xi to x'l). The number of chords does not depend 
on the choice of the maximal tree and is called chord number (or cyclomatic number). 
Whenever we add a chord to a maximal spanning tree we get a cycle, which is called 
fundamental cycle q. These chords are important because for T — )■ oo a general current 
of the form ([3]) can be written as a linear composition of the elementary currents trough 
the chords (see formula 0231) below). As an example, in Fig. [T]we consider a 4— states 
system, labeled by a, 6, c, where there are no transitions between b and c. We show a 
spanning tree where the two links left out are the chords (a, b) and (a, c). Moreover, the 
two fundamental cycles associated with each chord are displayed in Fig. [2J Note that 
this network has eight different spanning trees. 

Moreover, any cycle in can be written as a linear combination of the fundamental 
cycles, i.e., the fundamental cycles form an orthogonal basis (this was originally found 
by Kirchhoff [2H])- In order to see that we define the scalar product between a cycle 
and a pair of states as: 

{0 if (x, x') is not part of c 
1 if (x, x') is part of c and in the same direction . (18) 
— 1 if (x, x') is part of c and in the opposite direction 

This quantity gives an answer to the following question: is (x, x') part of the cycle c. 
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d = a 




d - a 




Figure 2. Graphical representation of relations (j20p and ((2T|) . On the left side we 
have the cycle C3 = (a, b, d, c, a), the first cycle on the right side is ci = (a, 6, d, a), and 
the second is C2 = (a, c, d, a). 



and if yes, are they oriented in the same direction? It is possible to show (see [21]) that 
any cycle c can be written as 



c = ^(c, {xi,x'i))ci. 



(19) 



The basic reasons for this formula are that a chord is not shared by two fundamental 
cycles and a chord must belong to a cycle c G G. This means that the increment of a 
cycle is given by 



i^(c) = ^(c, {xi,x',))K{ci). 



(20) 



Furthermore, due to the fact that In ^^||y^ is a particular case of a current K{c), as a 
corollary we get 



W{c) 



W(c) \>-\W(ci)J 



(21) 



The elementary current trough an edge (a:, x') can be written as a linear combination 
of elementary currents trough the chords I. Explicitly, from the finite time Kirchhoff 's 
law (|T5i) we obtain the asymptotic relation 



jT{x,x') = ^jr(xi,xO {ci,{x,x')) +0(^). 



(22) 



Therefore, from f[T^ we see that the general current (1131) takes the asymptotic form 
jT = l^Yl x'l) {ci, (x, x')) + 0{h 



\^Mxi,x'i)^f{x,x') {ci,{x,x')) + 0{^) 

I x,x' 

Y.Jt{x,,x[)K{ci) + 0{^). 



(23) 
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In particular, for the case of the entropy we have 



St = ^jTixi,x'i)\n 



W{ci) 



+ 0{-). 



(24) 



W(ci) 



Relation ( 12 2 p is normally written only for the (average) stationary current (see 
[TTtlSO]). however, we stress that it is valid for any stochastic trajectory. An elegant and 
heuristic proof of fl22|) for the stationary current can be found in [31], while a rigorous 
proof is given in [32] • 

A graphical representation of the decomposition of a cycle c into fundamental cycles 
using the scalar product (|T8|) for a 4— states system is showed in Fig. [2j The concepts 
developed here are central in studying the characteristic polynomial f ll2p . We proceed 
by deriving a general necessary condition for the symmetry to set in for a non-entropic 
current. As we show below, this condition is related to degeneracy in the set of values 
of the increments of cycles K{c). 

4. Necessary condition for the symmetry 

Using the Leibniz formula for determinants, the characteristic polynomial ( IT2|) can be 
written as 



where tt is a permutation of the states x G F and S'(F) is the permutation group 
associated with F. 

Moreover, each term of the determinant ( 125|) can be represented by a graph 
associated with the permutation vr on the network of states. Due to the bijectivity 
of TT, the graph is such that there is one transition rate entering and one transition 
rate leaving each state of the network. Therefore, as is well known in discrete 
mathematics [33], each permutation tt can be decomposed as product of disjoint cycles 
(with no common state). For example, in the case of seven states {1, 2, 3, 4, 5, 6, 7}, the 
permutation 7r{i) = 2i mod7, which is represented by {2, 4, 6, 1, 3, 5, 7}, can be written 
as (1,2,4, 1) o (3,6,5,3) o (7). We denote by C-,, the subset of the set of all cycles 
obtained from the cycles decomposition associated with tt. For example, in the case of 
TT{i) = 2i mod7 we have C,, = {(1,2,4, 1) , (3,6,5,3)}. Note that, as a consequence of 
the disjoint character of the cycles, the order of the cycles in is irrelevant. 

The increment related to C^r is written as 



(25) 



7re5(r) X 



cec^ 



(26) 
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and the transition rate is given by 

WiC^) = n Wic). (27) 

With these quantities, the characteristic polynomial fl2^ takes the expression 

Piz,y)= J2 sgn{n)exp{zK{C^))W{C^)g^{y) (28) 

7r65(r) 

where Qniy) is the term which comes from the fixed states of the permutation (vr(x) = x) 
and from the transpositions (7r(x) = x',7r{x') = x). Fixed states under the permutation 
contribute with the term —{y + A(x)), while transpositions contribute with the term 
Wx^x'Wx'^x- It then follows that 

9niy) = Y\ i-y-K^)) IJ (Wx^x'Wx'^x) ■ C^^) 

xer/'iT{x)=x {x,x'}e(r,r)/TT(x)=x' ,n{x')=x 

We define the set A = {C^,7r G 5'(r)} of the different disjoint cycle decomposition 
related to >S'(r). Equation fl28l) can be rewritten as 

P{z, y) = J2 i^KiC)) WiC) Yl sgn(7r)^^(|/), (30) 
ceA 7re5(r)/c^=c 

where the last sum is restricted to permutations vr such that = C. Since C^-i — C*7r) 
where C,, is the set of cycles composed of cycles in the reversed direction, we obtain 

P{z,y) = J^exp {zK(C)) W(C) Yl sgn(7r)^.(y) 

= 5^exp(-^i^(C))iy(C) J2 ^M^)9Ay) 

CeA 7r65(r)/c^_i=c 

= ^exp(-^K(C))H^(C) Y sgn(7r)(7.(y), (31) 
CeA 7r65(r)/c^=c 

where we used the properties K{C) = —K{C), (7^-1 (x) = g-jrix) and sgn(7r~-'^) = sgn(7r). 

Imposing the symmetry 

P{-E - z,y) = P{z,y) for all y, (32) 

and using formulas (15U]) and ([SI]), we get 

Y,^^PiEK{C))exp{zK{C))W(C) Yl ^M^)9Ay) (33) 
CeA 7r6S{r)/c^=c 

= Y i^^iC)) W{C) Y sgn(7r)^?.(y). (34) 

CeA nes{T)/c^=c 

We immediately see that, due to the independence of the family of exponential functions, 
the different and non- vanishing values of the set K{A) = {K{C),C G A} play an 
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important role in the analysis of this equation. In particular, we consider {K{ci)}: the 
subset of K{A) composed of cycle decompositions that contain just one fundamental 
cycle. If the following three conditions are fulfilled: 

• K{ci) for all fundamental cycles; 

• there is no degeneracy in {K{ci)}] 

• the set K{A) — {Ji (q)} is disjoint with{i^'(Q)}. 

Then relation (l34l) implies that 

exp {EK{ci)) = for all chords /. (35) 

In this case, we see from fl2^ and f l2^ that Jt = + 0{j;): we have the familiar 
GCEM symmetry because at large times Jt is just the entropy current (up to a constant). 
Therefore we arrive at the following conclusion: in order to have a symmetric current 
asymptotically different from entropy either at least one of the increments in the set 
{K[ci)} has to be zero or we need degeneracy in the increments of the fundamental 
cycles, in the sense that at least one of the last two conditions have to be broken. 
This is a main result of the present paper and we think it should play an important 
role in developing a more general theory for the symmetries of large deviation function 
associated with time-integrated currents. 

Two remarks are in time. First we stress that this is a necessary condition on 
the full spectrum of the modified generator: it could be that there is a situation 
where a non-entropic current for which this condition is not fulfilled (implying in an 
asymmetric characteristic polynomial) still has a symmetric maximum eigenvalue (which 
gives the scaled cumulant generating function). The other remark is a comparison of 
our work with [18] , where the condition f l35|) was obtained as a necessary and sufficient 
condition for the GCEM symmetry through a similar analysis of the characteristic 
polynomial ( 125|) . The "single macroscopic current" considered in [18] is proportional 
to the entropy current and here we are interested in non-entropic currents, which might 
have a symmetric characteristic polynomial precisely when condition (1351) is not satisfied. 
Another difference is that in [18] a sort of coarse-grained current is considered, where 
the number of cycles, and not of the elementary jumps, is counted. In this case one 
takes a derived Markov process which is obtained from the original process by cutting 
off cycles 

Even though we found a general necessary condition, this is still far from a sufficient 
condition. As we show next, in order for a symmetry to set in the transition rates have 
to fulfill some constraints, which depend on the kind of degeneracy we have among the 
increments in the set K{A). In the next section we perform a full analysis of all possible 
symmetric time-integrated currents different from entropy that arise in a 4— states 
system. Specifically, we consider the network of states of Fig. [H where all but one pair 
of states are connected. We show that in this case, two classes of symmetric currents 
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different from entropy exist, one of them corresponding to the symmetry previously 
found in [22] • Moreover, in the Appendix B we perform a similar analysis for the fully 
connected network, where there are several different classes of symmetric time-integrated 
currents. 



5. Three cycles and 4— states system 

The network of Fig. [1] has two 3— jumps cycles and one 4— jumps cycle. They are 
Ci = {a,b,d,a), C2 = {a,c,d,a), and C3 = {a,b,d,c,a), where C3 is the 4— jumps cycle. 
Besides the forward cycles there are also the three backward cycles, which are denoted by 
Ci, C2, and C3. The increment of the cycles is denoted by Ki and the forward (backward) 
rate by Wi (Wi), with i = 1,2,3. Moreover, for later convenience we write the escape 
rates from states b and c as A(c) = Ai and A(6) = A2. 

If we choose the spanning tree showed in Fig. [T]the two chords are (a, b) and (a, c), 
whereas the fundamental cycles are ci and C2, respectively. The rate and increment of 
the cycle C3 fulfills the following relations. From ( 120|) we obtain (see Fig. |2]) 

K3 = K,- K2, (36) 

and from fl2T|) we have 

W3 W1W2 



W3 W1W2 

Furthermore, from ( I23p . asymptotically, a general current is written as 



(37) 



jT = K,jT{a, b) + K2jT{a, c) + 0(i). (38) 
For the case of the entropy the above formula becomes 

st = kA In =Lj^(a, b) + In =^j^(a, c) + O(-) 39 
\ Wi W2 J T 

where is a constant. Hitherto we have defined entropy as the current for which 
Kg = this section we denote st any current that is asymptotically proportional to 
entropy. Note that in this case the GCEM symmetry (ITU]) changes to is{,z) = Is{—Es—z), 
where Es = l/K., (see [I9]). 

With the above definitions we can now look for the currents with a symmetric 
characteristic polynomial. As we showed in the previous section, the symmetry implies 
the fulfillment of equation f lM]) . In the present case this equation takes the form 

-(e^3(.+i^)]y3 + e-i^3(.+i^)|y3) 

= (e^^ + e-^i^WT) (-Ai -y) + (e^^W2 + e'^^^W^) (-A2 - y) 

- (e^^W3 + e-^^'^Ws) , (40) 
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where we used the fact that the permutations related to ci and C2 have positive sign 
and the permutation related to C3 has negative sign. We also showed in the previous 
section that in order to have a symmetric current different from entropy either one of the 
increments of the fundamental cycles has to be zero, or one of the fundamental cycles 
must have the same increment as another cycle (fundamental or non-fundamental). In 
the present case this leaves us with six possibilities: 

• = K2; 

• K, = 

• = or Ki = 0; 

• Ki = 2K2, which gives = K2, or K2 = Si^i, which gives = —Ki. 

We first consider the family of currents which have the cycle increments 

Ki = K2 = Ka, implying = 0. This means that all currents is this family have 
the following common asymptotic behavior, 

A"'^ = K^iMa, b) + jr(a, c)) + 0(i). (41) 

Note that the sufficient relation for the symmetry fHU]) depends only on the values of 
Ki , K2 and K^. Therefore, it is identical for all currents in this family. 

We want to find out under which conditions a current with the assymptotic form 
(HTj) is symmetric and different from st- From equation fl40p we obtain 

Wi + W2 



where we are using for the symmetric factor related to current fHTl) . This relation 
defines the value of Ea as a function of the transition rates. Moreover, also from (l40l) 
we get 

Wi + W2 WiXi + W2X2 



Wi + W2 WiXi + W2\2 



(43) 



This equation has two solutions: Ai = A2 and = = =. Clearly, from equations (139|) 
and ( 14T|) . if the second solution is satisfied we have = st/ Ea + 0{^) (with Kg = 1). 
Therefore, we obtain that the spectrum related to is symmetric and j!^^ is different 
from entropy if 

W^ Wo 

Ai = A2 and = ^ =■ (44) 

Wi W2 



Therefore, we have the symmetry Ia{z) = Ia{—Ea — z) with Ea given by (l42l) . 

The second family of currents j'^'^ corresponds to the case Ki = —K2 = K 
= 2Ki3. The asymptotic behavior of the currents in this family is 



= K^iMa, b) - Ma, c)) + 0(i). (45) 
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Following the same procedure of the previous case, from equation fl40l) we obtain that 
the symmetric factor is defined as 

eMKpEp) = \ (46) 

and the transition rates have to fulfill the constraint 
Wi + W2 WiXi + W2A2 



Wi + W2 WiXi + W2X2 



(47) 



This equation has the solutions Ai = A2 and = = =, where if second one is fulfilled j^'' 
becomes proportional to entropy. In addition, since is different from zero, relation 



40|) gives us an extra constraint, which is 

fWi + W2\^ W1W2 



— ^ ' M =^11:1^^ (48) 

\Wi + W2j W1W2 

where we used exp(2i^/3_E') = and relation (137|) . This equation has the solution 

WiWi = W2W2 and ^ = ^ (which gives entropy). Therefore, we conclude that Ifi^z) 
is symmetric, with the symmetric factor Ep given by fH6l) . and different from entropy if 

^^Z}^ - and (49) 

WiWi = W2W2 Wi W2 

For the case K2 = 0, we take = Ki = K . Hence, the associated family of 
currents have the asymptotic form KjT{a,b). Equation (HOj) gives the constraint 

Wi WiXi + W3 



Wi WiXi + W3 



(50) 



This is satisfied only if ^ = which gives =- = 1. This leads to the entropy current 
being proportional to the elementary current trough (a, b) and, therefore, there is no 
new symmetry in this case. The case i^i = is analogous. 

Finally, for the case Ki = 2K2 = K3 = K, equation fHOj) gives the constraint 

W2 W2X2 + W3 



W2 W2X2 + W3 

2 



(51) 



This is satisfied only if = which leads to = j . This last relation implies 
the associated current being proportional to entropy: there is no new symmetry. The 
case K2 = 2Ki is analogous. 



5.1. Example: two sites SEP 



In order to illustrate these symmetries we consider the symmetric exclusion process 
(SEP) [3H437] . This is a one dimensional transport model which is driven out of 
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equilibrium by the boundary dynamics. Particles enter and leave the system at the 
left and right boundary, while they diffuse in the bulk. They also interact in the bulk 
by imposing that the maximum number of particles per site is one. More generally, we 
will consider a SEP where the transition rates for particles to enter or leave the system 
depends on the bulk density. We will restrict our analysis to a two sites system which 
has the 4— states network showed in Fig. [TJ The transition rates for this two sites SEP 
are shown in Fig. [31 Note that, the subscript 1 (2) is related to transitions at the 
boundary with the other site being occupied (empty). 

Considering Fig. [1], we identify the state where the left (right) site is empty and 
the right (left) site is occupied with a (d) and the state where both sites are occupied 
(empty) with b (c). In the basis {a,b,c,d}, the generator ([2]) for the two sites SEP 
becomes 

/ — 1 — ai — ^2 cti So 1 \ 



ai - 


- S2 Ctl S2 


1 




-f3i - 5i 


Si 


72 


-72 - a2 


Ci2 


1 


7l 1^2 


-1-/32 



(52) 



\ 1 7l P2 -1 - /52 - 7i / 

One physical current that we are interested in is the current between the system and 
the left reservoir j^, which is asymptotically equal to minus the current between the 
system and the right reservoir. This is a functional of the stochastic trajectory such that 
whenever a particle enters (leaves) the system in the left boundary the current increases 
(decreases) by one. Therefore, for the two-sites SEP, the left boundary current is such 
that the only non-zero increments are f{a,b) = 1 and f{c,d) = 1. Hence, in the long 
time limit - where jr(c, d) = jTia, c) + 0{j;) and jrib, d) = jria, b) + O(^) - we have 

jlr = Ma, b) + jt(c, d) = Ma, b) + jria, c) + 0(i). (53) 

This current belongs to the family j^"^ with = 1. With the generator (I52p we have 
Wi = aiSi, Wi = 7i/3i, W2 = €(262, and W2 = 12^2- Therefore, the entropy fl39|l . with 
Kg = 1, reads 

ST = In ^Ma, h) + In ^jT(a, c) + 0{h. (54) 
7i/3i I2P2 J- 

An important point is that in general st and j'^ are not proportional, however for the 
standard SEP, that is, without dependence of the boundary transitions on the bulk 
density, the transition rates are such that and st are proportional [38]. 

Furthermore, for the two sites SEP the current from the left reservoir can be divided 
into two parts: the current when the right site is occupied, which is jT{a,b), and the 
current when the right side is empty, which is jT{c,d). We now consider a second 
current which is the difference between these two contributions: the current from the 
left reservoir when the right site is occupied minus the current from the left reservoir 
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Figure 3. Transiton rates for the two sites SEP with the boundary rates depending 
on the bulk density. 



when the right site is empty. More clearly, in the large deviation regime this current 
reads 



■Id 
Jt 



jria, b) - jric, d) = jria, b) - jria, c) + O(-), 



which is of the form j!^^ with K/^ = 1. 

The modified generators for these currents are obtained from fITT]) . As an example, 
we write down the modified generator for the current ( 155|) . which is 

^ — 1 — ^2 — a2exp{z) di 1 \ 

I32exp{z) -[52-^2 82 

7i — 7i — «! aiexp(— 2;) 

\ 1 72 /3iexp(z) -l-/3i-72 j 



(55) 



(56) 



Our theory predicts that the maximum eigenvalue of the modified generators 
associated with currents ( 153|) and (!55l) are symmetric, with symmetric factor given by 

and fl46p . respectively, if some constraints on the transition rates are satisfied. More 
clearly, if we consider 02 and 82 as depending on the other transition rates, from (jUj) 
we see that has a symmetric and non-entropic scaled cumulant generating function 
if 



a2 = 5i + (3i- 72 



and 



(57) 



Moreover from 



for if: the conditions for the symmetry are 



«2 = (^1 + /3i - 72 

^ '-012P212 



and 



ai5i _^ 02^2 



liPi I2P2 



(5^ 



In Fig. Hlwe show the scaled cumulant generating function, obtained by numerically 
calculating the maximum eigenvalue of the modified generator for the case ai = 1, 
5i = 0.31, f3i = 0.5, 7i = 0.1, 72 = 0.8, and /32 = 0.13. The scaled cumulant generating 
functions for the non-entropic currents are appropriately rescaled by their symmetric 
factors so that they touch the horizontal axis at —1. 
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Figure 4. Three symmetric currents: Ia{zEa) in black, Ip{zEfs) in red, and Is{z) in 
blue. We are considering the currents we defined for the two-site SEP (see text) which 
corresponds to Ka = Kp = = 1. The value of the transition rates are ai = 1, 
5i ~ 0.31, Pi ~ 0.5, 7i ~ 0.1, 72 ~ 0.8; and /32 — 0.13. Moreover, a2 and 82 arc given 
by relation (|57)) . 



Lastly, a very relevant observation is the following. For the present 4— states 
model the space of all possible asymptotic currents in the large deviation regime is two- 
dimensional: there are two independent elementary currents because the chord number 
is two. This space of currents can be represented in the Ki x K2 plane. In this plane 
the entropic currents, i.e., the currents that satisfies the GCEM symmetry, are in a line 
given by 



Our results shows that there are two more lines in this plane with symmetric currents: 
for the family of currents j^"^ we have K2 = Ki and for j^'^ we have K2 = —Ki. 
Whereas the slope of the line for the GCEM symmetry depends on the transition rates 
for the other two new symmetries this is not the case. Moreover, different from the 
GCEM symmetry, the a and /3 symmetries set in only when the transition rates fulfill 
some constraints. This situation is depicted in Fig. O 

6. Time-reversal as the origin of the symmetries 

As is well known, the GCEM symmetry is a direct consequence of the fact that entropy 
is the logarithm of the weight of the stochastic path divided by the weight of the reversed 
path. We argue that also the new symmetries we found in this paper have their origins in 
time-reversal. However, for non-entropic currents, this becomes apparent only when we 



K2 




(59) 
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Figure 5. The three symmetric currents m the f{a,c) x f{a,b) plane. The value of 
the transtion rates, which determine the slope of the line for the GCEM symmetry, 
are ai = 1, <5i = 0.31, /3i = 0.5, 71 = 0.1, 72 = 0.8, and = 0.13. 

group the trajectories in a certain way. This group of trajectories depends on the current 
that we are considering and it is determined by the degeneracies of the increments of the 
cycles. Specifically, we consider the 4— states system of the previous section and provide 
a handwaving demonstration of the symmetries of the family of currents j^^"^ and 
The presentation here is closely related to reference [22], where the symmetry for the 
current j!^^ was considered. The new results here are the symmetry for the current 
which comes from a different grouping, and the conjecture that the group of trajectories 
is determined by the degeneracies in the increment of the cycles. We remark that a 
(completely different) grouping of trajectories was also considered in [39], for the study 
of the effects of coarse-graining on the fluctuation relations. 

Before going into the grouping of paths, let us first present a simple demonstration 
of the fluctuation theorem. The weight of a stochastic trajectory X[o,r] is given by 

Af-l 

P[X[o,T]] = exp[-A(x(tA/))(T - ^m)] JJ exp[-A(x(ti))(ti+i - ti)]. (60) 

i=0 

The exponential of the waiting times multiplied by the escape rates come from the 
fact that we are dealing with a continuous time process and we should also multiply 
this expression by P(x(to),^o) which we are assuming to be uniform for simplicity. 
The reversed trajectory, where the system starts at state x(tAf) (also with a uniform 
probability distribution of sates) at time and jumps from state x(tj+i) to state x(ti) 
at time T — tj+i, is denoted by X^q^t]- Hence, 

M-l 

P[X[o,r]] = exp[-A(x(tM))(T - ^a/)] JJ ii'xCt.+o^xCti) exp[-X{x{ti)){ti+i - ti)]. (61) 

i=0 
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The entropy current ([5]) is related to the weight of a trajectory divided by the weight of 
the time-reversed trajectory by 

exp(^T[X[o,T]]) = ^^4- (62) 

This is a fundamental relation and it is at the origin of the infinite time fluctuation 
theorem we consider here as well as all other fluctuation relations. From fl62l) it follows 
that 



P(X[o,T]) exp{-A)6{ST[X[o,T]] - A) = F[X^o,t]]S{St[X[o^t]] - A) 

= nx^}]6{St[Xi^]] + A), (63) 

where we used the relation S'r[X[o,r]] = — 'S't[X[o,T]]- Summing the above relation over 
all trajectories and taking relation (|4]) into account we obtain 

In the limit T — oo the above formula implies the GCEM symmetry ([7]). A similar 
kind of demonstration also holds for currents of the form and This happens 
because a relation analogous to (162|) is valid for these two families. However, in order 
to see it we have to consider a functional of an appropriate group of trajectories. The 
demonstrations that follow are rather heuristic and we plan to provide a precise proof 
of these symmetries trough time-reversal of a group of trajectories in future work. 

6.1. The current j!^"^ 

We now consider the current Jt[-^[o,t]]; which is defined by the increments 
f(a, b) = /(., c) = m d) = /(c, = i In '^■+^'>;"- 

2 {Wi + W2)Wd^a 

/(d,a) = ln^^. (65) 

Wa^d 

It is easy to check that this current is in the family j^"-* with = In (which 
gives from E„ = 1). 

Let us now define the group of trajectories. Two trajectories ^[o,r] and X'q belong 
to the same group if they have the same number of jumps taking place at the same times 
tj. Moreover, if x{ti) = a then x'(tj) = a, and if x{ti) = d then x'{ti) = d. On the other 
hand, if x(tj) = b or x{ti) = c then x'(tj) = 6 or x'{ti) = c. In this way, a trajectory 
that stays in states 6 or c during m of the M + 1 time intervals is part of a group with 
2"^ trajectories. The group of trajectories is denoted by {X[o,r]}a, where the subscript 
indicates that we are dealing with a current of the type Its weight is just the sum 
of the weights of all the trajectories in the group, that is, 

n{X[0,T]}a)= Yl nX[0,T]). (66) 

^[0,T]^{^lO,T]}oi 
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Moreover, from the increments fl65|) . we see that the current Jr[^[o,T]] is invariant within 
the group, meaning that it takes the same value for any X^q^t] G {X[Q X]}a- Hence, if we 
want to write this current as a functional of the group of paths {X[o.T]}a we can define 

■JT[{^[0,T]}a] = Jt[X[o,t]] with X[o,T] £ {X[o,T]}a- (67) 

This relation, added to (j4]) and (l66l) . leads to 

p(Jt = a)= Y1 n{Xio,T]}a)S{JT[{X^o,T]U-A). (68) 

Furthermore, the expression for the weight of a group of trajectories P({X[o,T]}a) can 
be written as 

A/-1 

P({^[o,T]}a) = n E ^-iu)^-iu+r) exp[-X{x{ti)){ti+, - 1,)], (69) 

where the sum J2x{ti)€{x{ti)}a '^^^^ '^^^ term for x(ti) = a,d and is over the states b 
and c if the paths in the group have x(ti) = b,c. Similarly we denote by P({X[o,t]}q) 
the weight of the group formed by the reversed trajectories. Then, the increments of 
In !!^!^'°'^'!" t, for the case x(ti) = b, c, can be written as 

exp[-\{b){ti+i - tj)] + w^(t^^^)^cWc^x{ti+i) exp[-A(c)(tj+i - tj)] 
Wx{u+i)^bWb^x{u_^) exp[-A(6)(ti+i - ti)] + w^x(t,_i)^cWc^a;(ti+i) exp [-A(c) (t^+i - ti)] ■ 

(70) 

In opposition to the ratio of single paths, the time- dependence in the exponential waiting 
time distribution does not cancel out in general. However, for the case A (6) = A(c), which 
is condition ( H4|) . the increment ( I70|) becomes a time- independent term which reads 

Therefore, when condition f H^ is fulfilled, with the choice of increments f l65|) . we expect 
that asymptotically 

.xp(j.iw„,.,}j) = 5paM + o(i). (72) 

This means that despite the fact of the current J being different from entropy, when 
paths are grouped appropriately a relation analogous to ( 1^ in this new coarse-grained 
space (of trajectories) is valid. Finally, summing over all possible groups of trajectories 
we obtain 

^^^i^^ = exp(-A), (73) 
P{Jt = A) ' ^ ' 

which proves the symmetry of the current J (and then of all the currents of the family 
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This demonstration has two important features. First, the grouping of paths is such 
that paths in the same group differ by cycles that have the same increment and that 
is the reason why the current under consideration remains invariant within the group, 
as in relation (1671) . The second feature is that the constraint on the transition rates 
(jH]) has the effect of making the ratio fITOl) time- independent. As we see next these two 
properties are also shared by the class j!^^ . 



6.2. The current j'^^ 

The grouping of trajectories is more complicated in this case and we have to use a 
different strategy to build an heuristic demonstration. We take the current Jj-fXp^T]], 
which is defined by the increments 

1 W1+W2 
f{a, b) = f{a, c) = f{b, d) = f{c, d) = - In 

2 W1 + W2 

/(d,a) = l. (74) 

This is a current of the from ( H5|l with Kjs = In ^^^^^ (implying in Ejs = 1). 

In order to roughly define the group {X[o,t]}/3 let us consider a three jumps sequence 
between times [tj_i,tj+2]. For two trajectories in the same group we assume that the 
first and last state in this three jumps term are the same in both trajectories. However, 
the two middle states can be different. In order to group cycle ci with C2 and C2 with 
Ci, we choose the following four jumps sequences to get grouped together: 

(a, b, d, a) with (a, d, c, a) 

(a, d, b, a) with (a, c, d, a) 

{d, b, a, d) with (d, a, c, d) 

{d,a,b,d) with {d,c,a,d). (75) 

For example, this means that if we look at this three jumps piece of two trajectories in 
the same group, for one trajectory we could have (a, b, d, a) and for the other (a, d, c, a). 
The reason for the more complicated grouping in relation to the previous case is that 
two paths in the same group differs not only by a state c instead of a state b, but also 
the order of the middle states in the three jumps sequence is inverted. 

Similarly to the previous case we expect that the current Jr[X[o_T]] is invariant 
within the group, i.e., 

JT[{X[o,T]}t3] = Jt[X[o,T]] withX[o,T] e {X[o,T]}/3- (76) 



We define P({X[o,t]}/3) as the sum of the weights in the group. Let us take a three jumps 
term in the sum In !!!!^'°'^'!^ t, such that, for example, we have either the cycle (a, fe, d, a) 
or the cycle (a, d, c, a) in the trajectories in the group. This piece of trajectory should 
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contribute with 

Wiexp{-X{b)Ati - X{d)At2) +W2exp{-X{d)Ati - A(c)At2) 
W2exp(-A(c)Ati - X{d)At2)+Wi exp(-A(rf)Ati - X{b)At2) 



(77) 



where Ati = tj+i — and At2 = — ^i+i- If the first condition in ( HUi) . which is 
A (6) = A(c), is satisfied the above term becomes 



H^i exp(-A(&) Ati - X{d)At2) + W2exp{-X{d)Ati - X{b)At^, 
W2exp{-X{b)Ati - X{d)At2) + Wiexp{-X{d)Ati - X{b)At^, 



(7J 



In order to get a time-independent term we need the second condition in f H9|) to be 
fulfilled. Explicitly, for l^iH^i = W2W2 the term ([78]) becomes Noting that when 
WiWi = W2W2, apart from f{a,d), the increments f l71|) are equal to | In ^ we claim 
that 

Jl^(|-^[0,T]|/3j J- 

After multiplying by the delta functional and summing over all groups of trajectories 
we get the final result, 

^fc^ = exp(-A). (80) 
P(Jt = A) 

Once again, the degeneracy of the increment of the cycles determines the group of 
trajectories and the constraints on the transition rates have the effect of making the 
increments in the sum In pUj^'"'"^^!^^ time-independent. We conjecture that whenever the 



full characteristic polynomial f l25|) associated with a non-entropic current is symmetric, 
the symmetry comes from the time-reversal of a group of trajectories and the two 
properties just cited should hold. 



7. Conclusions 

In this paper we have studied the symmetries of large deviation functions associated with 
non-entropic currents in pure jump processes. The most general result we obtained is the 
necessary condition for the appearance of a symmetric current different from entropy. 
This condition is related to degeneracies in the increments of the fundamental cycles 
of Schnakenberg network theory and we believe that it is a good starting point for 
the development of a general theory for the symmetries of large deviation functions 
associated with currents. 

As an example, we studied 4— states systems, where symmetric non-entropic 
currents were found. In this case we saw that these symmetries set in when the 
increments of cycles are degenerate and when the transition rates fulfill a set of 
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constraints, which depends on the currents under consideration. Moreover, with this 
example we learned an important lesson: the symmetries in non-entropic currents 
come from the time-reversal of a group of trajectories. The degeneracies in the cycles 
increments determine the group of trajectories, which is such that trajectories in the 
same group differ only by cycles with the same increment. Furthermore, the constraints 
on the transition rates have the effect of making the logarithm of the sum of weights 
of the paths in a group divided by the sum of weights of the the reversed paths time- 
independent. 

The demonstration of the new symmetries as the time-reversal of a group of 
trajectories we presented here is still handwaving and restricted to the 4— states and 
three cycles system. We plan to develop a precise proof in future work by considering 
a coarse grained space of trajectories where the (exponential of the) non-entropic 
symmetric current becomes the ratio of the weight of the path and the weight of the time- 
reversed path. Another interesting direction for future work would be to find situations 
where the characteristic polynomial of the modified generator is not symmetric but the 
minimum eigenvalue is still symmetric. It could be that in such case, if it exists, the 
symmetry in the large deviation function is associated with the time-reversal of some 
most probable trajectory (or set of trajectories), which would dominate in a sum over 
all trajectories with a given constraint. 
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Appendix A. Generalization of Kirchhoff 's law 

Keeping the dependence on the stochastic trajectory X^q^t] explicit, equation (fT5l) reads 
5^jt(x,x')[X[o,t]] = 0(1). (A.l) 

We would like to make two points about this relation. 

• First, and less generally, for the typical behavior, we get the usual Kirchhoff's law 
for the steady state current. By ergodicity, the typical behavior of jT{x,x')[X[o^T]] 
is given by 

jrix, x')[X[o^T]] Jinv{x, X') = Pinv{x)w^^^> - pinv{x')w^>^^, (A. 2) 

where pinv and Jinv are the mean density and current in the invariant steady state. 
In this case f lA.l|) implies 

J2Jinv{x,x')=0, (A.3) 

x'^x 
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which is the usual Kirchhoff 's law. 

• More generally, in the (non typical) large deviation regime, the finite time 
Kirchhoff 's law flA.l|) imply that the large deviation function of the joint elementary 
current is infinite if the usual Kirchhoff 's law is not fulfilled [22] • That means, if 
we define the probability distribution of the family {jT{x,x'), {x,x') G F x F} as 

Pm = {a}) = J2^{X[^,T])llHMx,x') - a^^O , (A.4) 

Xt x,x' 

we have that the large deviation function associated with it respects the relation 
Hi'^}) = 00 if 3x such that a^x' 7^ 0. (A. 5) 

x'j^x 

In other words, the distribution of elementary currents around the atypical limit 
which do not fulfill Kirchhoff's law decreases faster than exponentially (and is not 
considered in this article). 



Appendix B. Fully connected 4— states network 



We now turn to the fully connected 4— states network, which is shown in Fig. IBli The 
states are now denoted be A, B, C and D. The set contains 14 cycles, which are: 

• four cycles with three states: Ci = (A, B, D, A), C2 = {A, D, C, A),C3 = {B, C, D, B), 
and C4 = {A,C,B,A); 

• three cycles with four states: C5 = {A, B,C, D, A),cq = {A, B, D,C, A),C'r = 
{A,D,B,C,Ay, 

• and the respective reversed cycles. 

Taking the spanning tree of Fig IBll the three chords are {A,B), {C,A) and {B,C). 
Whereas the three fundamental cycles are ci, C2 and C3, respectively. The increments 
of the four 3— states cycles are denoted by Ki, with i = 1,2,3,4. On the other hand, 
the increments of the three 4— states cycles are denoted by L,, with z = 5, 6, 7. The 
forward (backward) rates are represented by Wi (Wi), for i = 1,2, 3, 4, 5, 6, 7. Moreover, 
we define the escape rates as 

Ai = X{C), A2 = X{B), A3 = A(A), A4 = X{D). (B.l) 

We now have three independent elementary currents and, from equation (1251) . a 
general current has the asymptotic form 

Jt = K,jt{A,B) + K,jT{C,A) + K,jTiB,C) + 0{^), (B.2) 

which implies 

/ 1^1 Wo W-? \ 1 

st = kJ In =^MA, B) + In =^3t{C, A) + In =^3t{B, C) + O(-). (B.3) 

\ Wi W2 W2 J J- 
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Figure Bl. Fully connected network of states. On the right, the spannmg three is in 
blue and the chords in red. 



Furthermore, from relation (120 p we have 

''3 





= -K^ - K2 - 




-- Ki + 


Le = 


--K^ + K2 


Lj = 


-- K2 + K,, 


ED) J 


i;ives 




WJV2W3 




W1W2W3 




W1W3 




W1W3 


We 


W1W2 




W1W2 


W^ 


_ W2W3 


Wj 


W2W3 



(B.4) 



(B.5) 

Because in this case we have more cycles, there is a larger number of possibilities of 
grouping cycles with the same increment, and different groupings generate different 
restrictions in the transition rates: each specific grouping has to be treated separately. 
In the following we show two different cases where a symmetry different from the GCEM 
symmetry sets in. 

First we consider Ki = K2 = and K3 = M^. This gives a family of currents 
which have the asymptotic behavior 

1 , 



,(7) 
Jt 



K^MA, B) + AVt(C, A) + M,3t{B, C) + 0[ 



T' 



(B.6) 



Therefore, the set of increments becomes 

{i^i, i^2, i^3, ^4, L^. Le, L7} = {K^, K^, M^, -2K^ - M^, + M^, 2K^, + L^}, (B.7) 

meaning that also cycles C5 and C7 have the same increment. Depending on the values 
of and the grouping of cycles might change. We would like to consider the case 
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where only Ki = K2, = Lj and all other increments are different. In order for this 
to be true we need the following restrictions on and to be fulfilled, 

\K^\ ^ \M^\ 2\K^\ ^ \M^\ 

3K^ + M^^0 AK^ + M^^O + 2M^ ^ 2M^ + K^^0. (B.8) 

If all this constraints are satisfied, the equality (IMI) leads to the following constraints 
on the transition rates. Similar to the way we obtained equations (1421) and (1431) . if we 
consider cycles ci and C2 we get 

exp(K^E^) = = =^ = = = — . (B.9) 

^ ^ W1 + W2 W1X1 + W2X2 

The first equality defines the value of E^, while, as in equation ( H5]) . the solution of the 
second equality that leads to a current different from entropy is Ai = A2. Considering 
the cycle Cq we obtain 

--p'^^'^-'''-') -wwr [w;twJ • 

where the second equality comes from the first equality in (lB.9p . One solution for the 
above equation is W1W2 = W1W2, but in this case becomes entropy. The other 
solution, which is the one we are interested in, is WiWi = 14^2 W^2- Moreover, the cycle 
C3 gives 

eMM.E) = ^, (B.ll) 

which, with relations (IB.Sp . (1B.9P and the constraint WiWi = W2W2 gives 

J? = K, (^MA,B)+jT{C,A)+(\n^^ ' In C) j + 0(1). (B.12) 

Finally, the equations obtained from the increments of the cycle C4 does not bring any 
new constraints while the cycles C5 and cy give 

+ Wv 

exp(K^E^ + M^E^ = _ _ . (B.13) 

This equation leads to the additional constraint =-=- = =^^=-. Hence, we conclude 

that the (non-entropic) family of currents |jt^^| has a symmetric large deviation 
function, with symmetric factor E^ given by (IB.Op . if the following constraints on the 
transition rates are respected: 

\\ = A2 — 

WiWi = W2W2 and ^ _^ ^3_^4) 

m^W^ ^ Ws+We W2 

W2 W3 W5+W6 
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Let us make two remarks. The first is that in the above example we grouped cycles 
Ci and C2- However we could have chosen any pair from the four 3— states cycles, where 
each pair will give a different family of currents. Therefore, we have six different families: 
the one we treated and five more. In order to obtain the constraints for the other five 
one just have to follow the same procedure as above. The second remark is that if we 
break one of the conditions in fIB.Sp the degeneracies (and also the restrictions in the 
transition rates) change. For example, if we consider the case = 2L^ we have a set 
of increments given by {K^, K^, 2K^, —4:K^, 3K^, 2K^, 3K^} and (1B.6P becomes 

J?) = K,{MA,B)+MC,A)+MB,C)) + 0{^). (B.15) 

Since now cycles C3 and C5 are grouped together the restrictions on the transition rates 
are different from (lB.14p . More precisely, it can be shown that for = 2K^, besides 
the restrictions flB.14|) . we also need 



W3 W1W2 

in order to have the symmetry. 

As a second case we consider Ki = —K2 = Ks and K3 = Ms. This leads to the set 
of increments {Ks, —Ks, Ms, Ks, Ks + Ms, 0, —Ks + Ms} and the currents of this family 
have the asymptotic behavior 

jif) = KsMA, B) - KsMC, A) + MsMB, C) (B.17) 

Similarly to the conditions flB.14p . in order not to get any extra degeneracies the 
following restrictions on the increments have to be satisfied, 

Ks^O Ms^O 

\Ks\ ^ \Ms\ 2\Ks\ ^ \Ms\ \Ks\ ^ 2\Ms\. (B.18) 

We shall not present all the calculations for this case: they are similar to other 
calculations in section [5] and in this appendix. In the case that all the above conditions 
are fulfilled we obtain the following. Currents of the form 



j^^' =Ks\jT{A,B)-jT{C,A)+[\n^^j In =£^^'^(5, C) j + O(-), (B.19) 
have a symmetric large deviation function, with symmetric factor given by 



(B.20) 
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if the following restrictions on the transition rates are fulfilled: 

Xi = X2 
A3 = A4 
' W{Wi = W2W2 

W2 W3 Wa+y/W^W^^ 

This fully connected case gives us the following intuition. As soon as the network 
of states has a more sophisticated topology, with a larger number of chords, many other 
symmetric currents can be found. Therefore, we think that the finding of all possible 
symmetric currents, by using the methods of Sec. O becomes extremely difficult for 
larger and more complicated networks. 



and 



m W2 

Wi W2' 



(B.21) 
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